Power of the power-laws: lessons from unification of small and large time scales for 

evolution 
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We develop a "umfied' model that describes both "micro" and "macro" evolutions within a single 
theoretical framework. The eco-system is described as a dynamic network; the population dynamics 
at each node of this network describes the "micro" -evolution over ecological time scales (i.e., birth, 
ageing and natural death of individual organisms) while the appearance of new nodes, the slow 
changes of the links and the disappearance of existing nodes accounts for the "macro" evolution 
over geological time scales (i.e., the origination, evolution and extinction of species). In contrast to 
several earlier claims in the literature, we observe strong deviations from power law in the regime 
of long life times where the statistics is, usually, poor. 
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The recent surge in the modeling of biological evolution 
and extinction of species, using the concepts and tech- 
niques of statistical physics, has been stimulated partly 
by the claims (see [Q, |, ^ Q for reviews) that the statisti- 
cal distributions of several quantities associated with the 
extinction of species follow power laws. However, almost 
all of these models focus only on the "macro" evolution 
(i.e., the evolution of species on geological time scales). 
Neither the birth, ageing and, eventually, the death of the 
individual organisms nor the detailed population dynam- 
ics make explicit appearance in these theoretical descrip- 
tions. On the other hand, in reality, a species becomes 
extinct when its entire population is wiped out. 

Therefore, we develop a "wnz/iecf' model of an eco- 
system that describes both "micro" and "macro" evolu- 
tions. The eco-system is described as a dynamic network. 
The "micro" -evolution over ecological time scales, i.e., 
birth, growth (ageing) and natural death of individual or- 
ganisms is described by the dynamics within each node 
The network itself evolves slowly with time; over 
sufficiently long time scales populations of some species 
would drop to zero, indicating their extinction, and the 
corresponding nodes would be deleted from the network. 
On the other hand, appearance of new nodes, together 
with their own population of individual organisms, sig- 
nals origination of new species. In addition, the links 
of the network also change slowly to capture the adap- 
tive evolution of the species by random mutations over 
geological time scales. 

The dynamic network: At any arbitrary instant of time 
t the model consists of N(t) species each of which may 
be represented by one of the A'' nodes of a dynamic net- 
work; the total number of nodes is not allowed to exceed 
Nmax- Our model allows N{t) to fluctuate with time 
over the range 1 < N{t) < Nmax- The population (i.e., 
the total number of organisms) of a given species, say, i 
at any arbitrary instant of time t is given by ni(t). The 
limited availablity of resources, other than food, in the 
eco-system imposes an upper limit rimax of the allowed 
population of each species. Thus, the total number of 
organisms n{t) at time t is given by n{t) = 'Y^f^^ ni{t). 
Both Nmax and rimax are time-independent in the model. 



The interactions: Prey-predator interactions are cap- 
tured through the matrix J. The influence of species 
i on species i is given by J^ ; in general, Jij ^ Jji. The 
only restriction we impose initially on the elements of J 
is that Jii = 0, i.e., none of the organisms preys on any 
other member of the same species. Since in all practi- 
cal situations, food webs specify only the sign of Jy we 
allow the off-diagonal elements of Jij to take only the 
values and —1. Thus, Jij = 1 = —Jji indicates that 
j is the prey and i is the predator. Similarly, the situ- 
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respectively, to competition and cooperation between the 
species i and j. We assign the values -1-1 or —1 to the 
off-diagonal elemts of J randomly with equal probability 
in the initial state of the eco-system However, our 
model can be easily generalized to take into account any 
other architecture of food webs ^ ^ . 

We now argue that the matrix J accounts not only for 
the mier-species interactions but also miro-species inter- 
actions. First of all, note that if Jij > 0, then the species 
j is a prey of the species i if, simultaneously, Jji < 
whereas the species j cooperates with i if, simultane- 
ously, Jij > and Jji > 0. Therefore, if Jij > 0, the 
quantity {Jij — Jji)/^ is unity if the species j is a prey 
of the species i, but it vanishes if the species i and j 
mutually cooperate. Similarly, if Jij < 0, the quantity 
— {Jij — Jji)/2 is unity if the species j is a predator of 
i, but it vanishes if the species i and j compete against 
each other. Now, consider the two sums 
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where the superscript ± on indicates that the sum 
is restricted to only the positive (negative) elements Jij . 
The sum is a measure of the total food currently 
available to the z-th species whereas —S~ is a measure 
of the total population of the i-th species that would be, 
at the same time, consumed as food by its predators. If 
the food available is less than the requirement, then some 
organisms of the species i will die of starvation, even if 
none of them is killed by any predator. This way the 



2 



matrix J can account for the shortfall in the food supply 
and the consequent competition among the organisms of 
the species i . 

The collective characteristics of species: The age of an 
arbitrary individual organism, say, a of the species i at 
time t is denoted by the symbol X{i, a; t). In our model 
each species i is collectively characterized by & |^ : 
(i) the minimum reproduction age X^epii), (ii) the birth 
rate M{i), (iii) the maximum possible age Xmax{i), and 
(iv) the elements and Jji {j = 1,2, N). An in- 
dividual of the i-th species can reproduce only after at- 
taining the age Xrepii)- Whenever an organism of i-th 
species gives birth to offsprings, M{i) of these are born si- 
multaneously. None of the individuals of the i-th. species 
can live longer than Xmax{i)- Thus, even if an individ- 
ual manages to escape its predators, it cannot live longer 
than Xmaxi"!-) because of "natural death" caused by age- 
ing. 

The dynamics of the eco-system: The state of the system 
is updated in discrete time steps as follows: 
Step /- Birth: Assuming, for the sake of simplicity, the 
reproductions to be asexual, each individual organism 
a {a ~ l,...,ni{t)) of the species i {i — 1,2, ...N{t)) 
is allowed to give birth to M{i;t) offsprings at every 
time step t with probability (per unit time) pb{i,a;t) 
which is non-zero only when the individual organism age 
X{i, a; t) > Xrepii; t). 

Step H- Natural death: At any arbitrary time step t 
the probability (per unit time) of "natural" death (due 
to ageing) of an individual organism a of species i is 
Pdii,a;t). 

Step III- Mutation: With probability Pmut per unit time, 
all the species simultaneously re-adjust one of the incom- 
ing interactions Jy by assigning it a new value of either 

or —1 with equal probability 
Step IV- Starvation death and killing by prey: If 5^ is 
larger than S~ then food shortage will be the dominant 
cause of premature death of a fraction of the existing 
population of the species i. On the other hand, if > 
Ui — , then a fraction of the existing population will be 
wiped out primarily by the predators. In order to capture 
these phenomena, at every time step t, in addition to the 
natural death due to ageing, a further reduction of the 
population by 
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is implemented where ni{t) is the population of the 
species i that survives after the natural death step above. 
C is a constant of proportionality that basically sets the 
time scale of the population dynamics. If implementation 
of these steps makes < 0, species i becomes extinct. 
Step V- Speciation: At each time step, the niches (nodes) 
left empty by extinction of species are re-filled by new 
species, with probability psp. All the simultaneously re- 
filled nodes of the network originate from one common 
ancestor which is picked up randomly from among the 
surviving species. All the interactions Jy and Jji of the 
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FIG. 1: The total number of species N{t) is plotted 
against time; the corresponding parameter set is: Nmax = 
100, n^a, 1000, C = 0.1, = 0.001, p™„t = 0.001. 



new species are identical to those of their common an- 
cestor; each new species, however, either competes or 
cooperates with its ancestor species. The characteristic 
parameters Xmax, Xrep, Mj of each of the new species 
differ randomly by ±1 from the corresponding parame- 
ters for their ancestor. 

Probability of birth: We chose the time- dependent proba- 
bilty of birth per unit time as 



Xmax{i) - X{i,a) 



Xmax{'^^ Xrepi^^ 
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Pb{i,a) = 

iff X{i,a) > Xrepii) and Xmax ii) > Xrepii) (3) 



Note that in the limit of vanishingly small population, 
i.e., Hi 0, we have pbii,a) ^ 1 if Xii,a) — Xrepii) 
and, then, pi, decreases linearly Q as the organism grows 
older. However, since the eco-system can support only a 
maximum of n^ax individual organisms of each species, 
Pbii, a;t) ^ Q as n^it) rimax, irrespective of the age of 
the individual organism a |jll[] . 

Probability of natural death: We assume the probability 
of "natural" death (due to ageing) to have the form |12] 
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Xii,a)Mil)- Xrepii) 



Xmaxii)Mii) - Xrepii) 
if Xii, oi) > Xrepii) 



Xrepii)Mii) ~ Xrepii) 



Xmaxii)Mii) - Xrepii) 

if Xii, a) < Xrepii) 
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provided Xmaxii)Mii) > Xrepii)- In all other situations, 
Pdii,a) = 1. Note that, for a given Xmax and Xrep, 
the larger is the M the higher is the pd for any age X. 
Therefore, each species have a tendency to increase M 
for giving birth to larger number of offsprings whereas 
the higher mortality for higher M opposes this tendency. 

The longest runs in our computer simulations were con- 
tinued upto a maximum of five million time steps. If each 
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FIG. 2: The total number of organisms n{t) is plotted against 
time over a relatively short part of the evolution of the eco- 
system. The parameter values are identical to those in fig^. 



time step in our model is assumed to correspond to a real 
time of the order of one year, then the time scale of 5 
million years, over which we have monitored our model 
eco-system, is comparable to real geological time scales. 

Since we faced difficulty in getting high quality data, 
with reasonably good statistics, for Nmax > 100 and 
nmax > 1000, we have carried out most of our simula- 
tions with Njnax = 50, 100 and Umax = 100, 1000 only. 
The data obtained from the different runs, each starting 
from a random initial condition, were averaged. Both 
CRAY-T3E and SUN workstations were used for the sim- 
ulations. 

In figjl we plot the total number of species, N{t), in 
a particu ar run, starting from a single initial condition, 
upto half a million time steps. In figj^ we plot the corre- 
sponding variation of the total population n{t) over rel- 
atively short interval of 20, 000 time steps only. These 
clearly demonstrate that the evolution has periods of 
"stasis" during which organisms populations keep fluc- 
tuating; the stasis are interrupted by occasional bursts 
of rapid extinctions followed by slower recovery. 

The average distributions of the lifetimes of the species 
are plotted in fig.|| for one set of values of the parame- 
ters. Clearly, the data are consistent with a power-law; 
the effective exponent, which is, approximately, 2, is also 
consistent with the corresponding estimate quoted in the 
literature ||]. However, in fig.H the power law holds 
only over a limited range Since real eco-systems 

are much more complex than our model eco-system and 
the available fossil data are quite sparse, it is question- 
able whether real extinctions follow power laws and, if 
so, over how many orders of magnitude. 

We have also observed (not shown) that the higher is 
the mutation probability Pmut the lower is the lifetime; 
this is consistent with the intuitive expectation that the 
higher rate of mutation leads to higher levels of biological 
activity in the eco-system thereby leading to the extinc- 
tion of larger number of species. But, Psp had weaker 
effect on the same data. However, if psp is too small 
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FIG. 3: Log-log plots of the distributions of the lifetimes of 
the species in an eco-system with (a) Umax = 100, (b) Umax = 
1000. In (b) the data beyond the lifetime of 10 are not shown 
to emphasize the initial power law regime (symbolized by the 
straight line with slope —2 . The other common parameters 
for both the figures are Nmax — 50, psp = 0.01, pmut = 0.001. 
The symbols x and * in (a) correspond to maximum 

simulation times 5 x 10^, 5 x 10"', 5 x 10^ and 5 x 10®, respec- 
tively while the symbols x and * in (b) correspond to the 
maximum simulation times 10*, 5 x 10* and 4 x 10^, respec- 
tively. Each of the data points has been obtained by averaging 
over 18 to 176 runs, each starting from a new random initial 
state. 



to maintain adequate pace of recovery of the eco-system 
after mass extinctions, the entire eco-system collapses. 

Fig.^ shows the time-averaged age-distribution in the 
populations of a species as well as the distributions of 
Xmax, Xrep and M. We see that the minimum age of re- 
production Xren IS quitc Small, as usual in the employed 
ageing model The age distribution decays stronger 
than a simple exponential, indicating a mortality increas- 
ing with age as it should be . The genetic death ages 
Xfnax — 72 are far above the upper end ~ 31 of the 
age distribution, as is appropriate for animals in the wild 
I p^ . Finally, fig.|| shows the distribution of M(i); this is 
relatively much broader than the distributions of X^ax 
and Xrep- 

In this letter we have presented a unified model which 
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FIG. 4: Semi-log plot of the distributions of X (+), Xmax 
(x) and Xrep (*); the parameter values are same as those in 
figfl 



describes the birth, ageing and death of individuals and 
population dynamics on short time scales as well as the 
long-time evolution of species. Not only the total number 
of species and the inter-species interactions but also the 
collective characteristics, namely, Xrep, X^ax and M, of 
each species vary following a stochastic dynamics. Thus, 
our model is capable of self-organization. To our knowl- 
edge, there are only a few earlier evolutionary models 
based on inter-species interactions which describe pop- 
ulation dynamics of each species. The population dy- 
namics within the framework of Lotka-Volterra equations 
have been considered earlier for only a few species. 
But, these do not account for the age distributions as 
the entire population of each species is represented col- 
lectively by a single dynamical variable in contrast to the 
explicit birth, growth and death of individual organisms 
captured in our model. 
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FIG. 5: The distribution of M; the parameter values are same 
as those in fig|l|. 



Since we have observed strong deviations from power- 
law in the distributions of lifetimes of species even in 
the Sole-Manrubia model ||l5|, in spite of the absence 
of detailed "micro" -dynamics in the latter, we strongly 
believe that this is a generic features of evolution and 
extinction of species. It would be interesting to investi- 
gate the geographical effects on our model eco-system by 
re-formulating it on a lattice in the same spirit in which 
some lattice models of prey-predator systems have been 
formulated ilql. 
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